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We show that the Abell Cluster A586 exhibits evidence of the interaction between dark matter 
and dark energy and argue that this interaction implies a violation of the Equivalence Principle. 
This violation is found in the context of two different models of dark energy-dark matter interaction. 
We also argue, based on the spherical symmetry of the Abell Cluster A586 that skewness is not the 
most general quantity to test the Equivalence Principle. 
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Introduction. It has become rather consensual that the 
problem of the nature of dark energy and dark matter 
(hereafter DE and DM, respectively) is crucial in con- 
temporary cosmology. Even though, observational data 
is fully consistent with the ACDM parametrization, in 
order to get a deeper insight into the nature of DE and 
DM one must consider more complex models and, in par- 
ticular, the interaction of those components. However, 
so far no observational evidence of this interaction has 
been presented. In this work, we argue that study of the 
Abell Cluster A586 exhibits evidence of the interaction 
between DE and DM. Furthermore, we show that this in- 
teraction implies a violation of the Equivalence Principle 
(EP) . Our results are obtained in the context of two dis- 
tinct phenomenologically viable models for the DE-DM 
interaction. We consider the generalized Chaplygin gas 
(GCG) model [2], a unified description of DE and DM, 
where interaction is an automatic feature of this proposal, 
but also a less constrained approach where DE and DM 
are regarded as two independent components, but inter- 
acting (see e.g. [§])• We show that interaction between 
DE and DM implies a violation of the EP between DM 
and baryons which can be detected in self-gravitating sys- 
tems in stationary equilibrium. For sure, the EP - that 
is, the universality of free fall - is one of the cornerstones 
of general relativity, however its validity at cosmological 
scales has never been directly tested (see [3| and refer- 
ences therein). The EP can be expressed in terms of the 
bias parameter, b, defined as ratio of baryon over DM 
density, at a typical clustering scale (Mpc). Should the 
EP hold, b would be a constant over time since then all 
clustering species would fall equivalently under the action 
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of gravity. Inversely, clustering should reflect the viola- 
tion of the EP through a different behaviour for both 
species. Interaction between DM and DE induces a time 
evolution of b. 

In this work we shall focus on the effect of interaction 
on clustering as revealed by the Layzer-Irvine equation. 
Given that the EP concerns the way matter falls in the 
gravitational field, considering the clustering of matter 
against the cosmic expansion and the interaction with 
DE seems to be a logical way to test its validity. In what 
follows we shall see that DE-DM interaction implies a 
departure from virial equilibrium. First, we will set the 
formalism to address the DE-DM interaction and con- 
sider two phenomenologically viable models: one based 
on a ad hoc DE-DM interaction Q, the other in the GCG 
with an explicit identification of the DE and DM com- 
ponents Our observational inferences are based on 
the Abell Cluster A586 given its stationarity, spherical 
symmetry and wealth of available observations [5||. 

Quintessence model with DE-DM interaction. The 
Bianchi identities with coupling £ give origin to the fol- 
lowing homogeneous energy conservation equations: 

pDAI + ZHpDM = QHpDM , (1) 



p DE + 3Hp DE (l + ude) = -QHp DM ■ (2) 

Notice that these equations imply that the energy ex- 
change between DE and DM is adiabatic (see e.g. @ 
and references therein). Moreover, the basic assumptions 
used in these equations are a constant equation of state 
parameter pbe = ^dePde and the following scaling with 
respect DM energy density 
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for a constant r\. From the time derivative of Eq. J3]) 
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inserted into Eq. ((2|) together with Eq. (TJ yields: 

. _ (?? + 3^_d_e)^.d.Eo 
Q,de + ^DM a^ 7 ' 

The solution of Eq. (TJ is given by 

a pDM e Jl a 

a~ 3 pDM ^P-DEoOP + Odm 



(4) 
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The DE evolution is then derived from the scaling di- 
rectly, or from Eq. ([2j) combined with the scaling: 
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In this model, from Eq. (f5]) one can see that the bias 
parameter depends on time as follows: 
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The GCG model. Let us now consider the GCG model 
with an explicit identification of DE and DM, as dis- 
cussed in Ref. The GCG model is considered here 
as it fares quite well when confronted with various phe- 
nomenological tests: high precision Cosmic Microwave 
Background Radiation data Q, supernova data si. fol. . 
gravitational lensing [ill], gamma-ray bursts [l2|, cosmic 
topology IIJ and time variation of the electromagnetic 
coupling [141 ]. In order to obtain a suitable structure for- 
mation behaviour at linear approximation, lode = —1 
(see [| and references therein). For the GCG admixture 
of DE and DM, the equation of state is given by [l|: 



A 



(8) 



where A and a are positive constants. From Q, the DM 
and DE expressions for a constant DE equation of state 
are given by 



PDM =PDM a 



PDE =PDE 



-3(l+a) 



n D Eo+^DM„a-^ 1+ ^ 



n DEo + fW a- 3(1+Q) 



(9) 
(10) 



so that we recover Eq. (fS]), but now with 77 = 3(1 + a) 
and lode = — 1. 

The Generalized Layzer-Irvine equation. Let us now 
turn to the cosmological gravitational collapse and its im- 
plication for the EP. The core of our approach lies on de- 
viation of the classical virial equilibrium, in its standard 
Layzer-Irvine equation form. We argue that A586 data 
allows to establish this departure independently of the 
DE-DM interaction model considered. It is possible to 
write the energy density conservation for non-relativistic 



self-gravitating dust-like particles through the so-called 
Layzer-Irvine equation (TJ)]. The kernel of the method 
is to consider the Newtonian kinetic energy, K, per unit 
mass, while keeping the average momentum and mass, 
M. constant: 
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where a is the scale factor of the Robertson- Walker met- 
ric. It then follows that: 



p K = MdK/dV = d{MK)/dV cx a 



(12) 



It is assumed that the mass evolution of the cluster re- 
mains constant over the course of the observation. The 
energy transfer between DM and DE is negligible at this 
point. The potential energy per unit mass derives from 
the definition of the auto-correlation function, [l5| 



(13) 



W = —2nGa pdm / dr£(r)r 



where we have replaced the background energy density 
by the DM energy density. After considering the DE-DM 
interaction, it follows that 



W OC a 2 + dln Po M /dlna 



and hence 



Pw = MdW/dV = d(MW)/dV cx a 



(14) 



(15) 



This is the source of difference from the usual dust case. 
The Layzer-Irvine equation for the energies per unit vol- 
ume is just a chain rule of time derivative for the energy 
density where the time is parameterized by the scale fac- 
tor, hence: 



d d 

-t: (pdm) = a— (pdm) = - [2pk + (1 - C)pw] H . 
at oa 



from which follows 



Pdm + (2pk + Pw) H = C,pwH , 



(16) 



(17) 



where H = a/a is the expansion rate. 

Furthermore, writing in terms of the virial equilibrium 
factor 2px + Pw and the departure to static equilibrium, 
due to the DE-DM interaction, Eq. lfl7|) becomes 



pDM + H(2p K +p w ) 



(77 + 3lj de )H 
1 + n D M /tt D E a~ rl 



Pw- (18) 



As before, it is possible to see from the equivalent of Eq. 
((H for the GCG model (for which co DE = -1 0) that one 
can map Eq. (JTSj) for the generic interaction model into 
the GCG model via the relationship 77 = 3(1 + a). Next 
we shall apply these equations to the stationary Abell 
Cluster A586 for which px and pw can be computed, 
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so as to compare with the observed local measurements 
with the homogeneous-spawned interaction term: 



2pK + Pw — (Pw ■ 



(19) 



The Abell Cluster A586. In order to estimate the cou- 
pling between DE and DM from Eq. (fl"9]) one has to find 
a particular cluster to compute px and pw • It is conve- 
nient that the cluster is as spherical as possible and close 
to stationary equilibrium. Under these conditions, one 
can approximate the kinetic and potential energy densi- 
ties as: 



d K 
Pk =M—K ~ M— 



9 M 



dV 



V 



Cluster 2 



^ ^Cluster 



, r (i T „ ,,W 3 G 

Pw =M—W ~ M— ~ - 



M 2 

m Cluster 



dV 



V 



8ir < R> R 3 cluster 



(20) 
(21) 



where Mciuster and Rciuster are the cluster's total mass 
(galaxies, DM and intra-cluster gas) and radius, a v is 
the velocity dispersion as determined globally from weak 
lensing, and < R > is the mean intergalactic distance 

The cluster must be also relaxed, since the core of our 
method consists in estimating the EP violation from a 
deviation from the standard form of the cosmic virial 
theorem defined by Eq. (fl7|) set with no interaction. 

Given these constraints a particularly suitable cluster 
for our purpose is the Abell cluster A586 It is found 
that the mass profile in this particular cluster is approx- 
imately spherical and that it is a relaxed cluster, since 
it has not undergone any important merging process in 
the last few Gyrs The agreement between dynamical 
(velocity dispersion and X-ray) and non-dynamical mass 
estimates (weak-lensing) indicates that A586 is in fact a 
relaxed cluster. 

Considering that gravitational weak lensing is indepen- 
dent from equilibrium assumptions about the dynamical 
state of the cluster, it turns out to be the best mass es- 
timator. Therefore, in our analysis we assume [jj: 



Cluster 



(4.3 ± 0.7) x 10 14 M, 



(22) 



which corresponds to the total mass inside a 422 kpc 
radius region estimated using gravitational weak lensing. 

In order to have a coherent set of data, we take for the 
velocity dispersion [jj: 



<j v = (1243 ± 58) kms~ 



(23) 



as computed from gravitational weak lensing measure- 
ments. 

The mean intergalactic distance is estimated using the 
coordinates (right ascension-a c and declination-5 c ) of the 
31 galaxy sample provided in Ref. []|. Given that 
weak gravitational lensing data concerns a 422 kpc ra- 
dius spherical region and the 31 galaxies lie within a 
570/iyQ 1 kpc region, one has to select from the original 
sample the galaxies that lie within the range of interest. 
Since at the cluster's distance, one arcsecond corresponds 



to 2.9 kpc, we select from our sample the galaxies that 
have a c and 5 C such that: 
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(24) 



where a center and S cente r are the coordinates of the cen- 
ter of the cluster and A max = 145" is the angular di- 
mension corresponding to a radius of 422 kpc. From 
this procedure, we build a sub-sample containing 25 
galaxies. From this sub-sample coordinates one can es- 
timate the mean intergalactic distance by elementary 
trigonometry, the distance between any two galaxies i 
and j with coordinates (a c ,, £ c ,) and {ae c j, <^cj) is given by 



r~ = 2d 2 [1 - cos(a c 



a C j)cosS C iCos5 C j — sinS C iSin8 c 



where d is the radial distance from the center of the clus- 
ter to Earth. Therefore the mean intergalactic distance 
< R > is 



< R >- 



N . 

N gal (N gal -V |^ 



gal I 



(25) 



where N ga i is the number of galaxies in the sample. In 
our sub-sample, N ga i = 25 and hence we get the estimate 
< R >= 309 kpc. Using Eqs. {22]), ([23]) and < R > we 
can estimate the kinetic and potential energy densities, 
Eqs. (2SJ) and {2TJ: 



(26) 
(27) 



p K =(2.14 ± 0.55) x 10~ 10 Jm~ 3 , 
p w =(-2.83 ±0.92) x l0- 10 Jm- 3 



where the errors were computed using linear error prop- 
agation. 

It is worth mentioning that 



PK 

Pw 



-0.76 ±0.05, 



(28) 



instead of —0.5 as one would expect for a relaxed cluster 
considering the standard form of the cosmic virial theo- 
rem and no DE-DM interaction. 

DE-DM interaction and putative evidence of violation 
of the EP. In what follows we use our estimates of Pk and 
pw, Ec[s. P6]l. (|27| . and the latest cosmological WMAP 
data [FT} to show the evidence of DE-DM interaction. We 
also demonstrate that this interaction implies a violation 
of the EP between DM and baryons. 

Let us first look at the quintessence model with DE- 
DM interaction. From Eqs. ^ and J2J the DE-DM in- 
teraction is exhibited through a non- vanishing £ or equiv- 
alent^, from Eq. (]4]), by the condition rj ^ ~3lode- 

Thus, assuming that lode = — 1, &de = 0.72, 
Oflji/j = 0.24, one can estimate rj for which Eq. lfl9]) 
is satisfied for the redshift of the A586, z = 0.1708. We 
find that: 



rj = 3.82 



+0.18 
0.17 



(29) 



Thus, since Eq. ((29]) satisfies the condition rj ^ 
— 3wd_e, one concludes that DE and DM are interact- 
ing. Notice that, as observations suggest a recent DE 
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dominance, then £ < 0, and from there follows that 
rj > — 3ujde- This means that Eq. lj29f not only suggests 
that DE and DM are interacting, but also, as expected, 
that the energy transfer flow is from DM to DE. 

Let us now turn to the CGC model. With the iden- 
tification of components suggested in DE-DM inter- 
action is expressed by the condition a^O. In order to 
see the effect of interaction in the GCG model, we pro- 
ceed as before using Eqs. (fl9|). ([26| and f2"7j) . from which 
follows: 

a = 0.27lg;gg. (30) 

Thus, the condition a ^ holds, meaning that the 
A586 data is consistent with the identification of compo- 
nents suggested in [3] for the CGC model. Notice that 
for a = the GCG model corresponds to the ACDM 
model. Moreover, it is interesting to point out that the 
value a ~ 0.27 is approximately consistent with values 
found to match the bias and its growth from the 2dF 
survey (see [3] and references therein). 

Evidence on a possible violation of the EP implies the 
time dependence of the bias parameter. We depict in 
Fig. [TJ the evolution with redshift of the normalized bias 
parameter predicted by Eq. (7]) where only gravitational 
effects are considered. Even though other astrophysical 
effects might affect the way DM and baryons fall under 
gravity, for EP purposes, gravity is the only relevant in- 
teraction that offers a clear drift on a cosmological time 
scale. Clearly, one expects that for large samples those 
effects would average out for non-cosmological drifts and 
thus lead to possible detection in large surveys. 
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Figure 1: Normalized gravitationally induced bias parameter 
as a function of the redshift, where 615 = b(z = 15), z = 15 
being a typical condensation time. 



Figure Q] shows that &(z)/&i5 has undergone a sharp 
change in the recent past, a clear signal of the violation 
of the EP due to the DE-DM interaction. This abrupt 
variation corresponds to the period when energy transfer 
from DM to DE becomes significant [z ~ 0.5). 

Discussion and Conclusions. In this work we have ar- 
gued that the properties of the A586 suggest evidence of 
the interaction between DE and DM. We stress that the 
considered models to describe the DE-DM interaction are 
consistent with known phenomenological constraints. We 
have also argued that this interaction does suggest a vio- 
lation of the EP that should be detectable in large scale 
cluster surveys via the analysis of the time dependence 
of the bias parameter. We find that this violation is in- 
dependent of the interaction model between DE and DM 
and entails a redshift evolution of bias parameter given 
by Eq. j7j and depicted in Figure [TJ Our conclusions are 
independent of the DE-DM interaction model, generic 
or GCG. Actually, a violation of the EP is reported to 
be found in other DE models [l8j]. For the GCG model 
we find that the detection of interaction precludes the 
ACDM model (a = 0). Furthermore, the obtained value 
for a is approximately consistent with results for the bias 
and its growth obtained by the 2dF survey Q. Consis- 
tency of our results with observational data concerning 
interaction and further implications of the detected 
interaction between DE and DM, for instance, in what 
concerns the motion of the satellite Sagittarius galaxy 
[2(ij |. are discussed in [2l| . 

It is interesting to point out that our results indicate 
evidence for violation of the EP between baryons and 
DM using data extracted from the A586, a notoriously 
relaxed and spherically symmetric structure. This seems 
to imply that the suggestion that cosmological evidence 
for a violation could be detected via skewness [l(| does 
not hold. Indeed, spherical symmetry implies that skew- 
ness vanishes given that it is an odd parity spatial func- 
tion. Thus, while the virial equilibrium may in principle 
reveal the violation of EP due to the DM-DE interaction, 
skewness is unable, by definition, to detect it in this par- 
ticular symmetry. The spherical symmetry of A586 and 
our detection of violation of the EP via virial equilibrium 
exemplifies this point. 
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